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NGC 1569 is a nearby dwarf irregular galaxy which underwent an intense 
burst of star formation 10 to 40 Myr ago. We present observations that reach 
surface brightnesses two to eighty times fainter than previous radio continuum ob- 
servations and the first radio continuum polarization observations of this galaxy 
at 20 cm, 13 cm, 6 cm, and 3 cm. These observations allow us to probe the rela- 
tionship of the magnetic field of NGC 1569 to the rest of its interstellar medium. 
We confirm the presence of an extended radio continuum halo at 20 cm and see 
for the first time the radio continuum feature associated with the western Ha arm 
at wavelengths shorter than 20 cm. Although, in general, the spectral indices de- 
rived for this galaxy steepen as one moves into the halo of the galaxy, there are 
filamentary regions of flat spectral indices extending to the edge of the galaxy. 
The spectral index trends in this galaxy support the theory that there is a con- 
vective wind at work in this galaxy. There is strong polarized emission at 3 cm 
and 6 cm and weak polarized emission at 20 cm and 13 cm. We estimate that 
the thermal fraction is 40-50% in the center of the galaxy and falls off rapidly 
with height above the disk. Using this estimate, we derive a total magnetic 
field strength of 38 /iG in the central regions and 10-15 /iG in the halo. The 
magnetic field is largely random in the center of the galaxy; the uniform field is 
~ 3-9 /iG and is strongest in the halo. Using our total magnetic field strength 
estimates and the results of previous observations of NGC 1569, we find that the 
magnetic pressure is the same order of magnitude but, in general, a factor of a 
few less than the other components of the interstellar medium in this galaxy. The 
uniform magnetic field in NGC 1569 is closely associated with the Ha bubbles 
and filaments. We suggest that a supernova-driven dynamo may be operating in 
this galaxy. Based on our pressure estimates and the morphology of the magnetic 
field, the outflow of hot gas from NGC 1569 is clearly shaping the magnetic field, 
but the magnetic field in turn may be aiding the outflow by channeling gas out of 
the disk of the galaxy. Dwarf galaxies with extended radio continuum halos like 
that of NGC 1569 may play an important role in magnetizing the intergalactic 
medium. 

Subject headings: galaxies: individual (NGC 1569) — galaxies: ISM — galaxies: 
irregular — galaxies: magnetic fields — galaxies: starburst — radio continuum: 
galaxies 
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Introduction 



The structure of a galaxy's interstellar medium (ISM) is shaped by its stars. Stars 
dissociate and ionize the surrounding cocoon of cold gas from which they formed. Winds from 
massive stars and supernovae explosions also create bubbles of hot gas. The combined effects 
of multiple stars in a cluster create larger bubbles of hot gas referred to superbubbles. The 
low mass of dwarf irregular galaxies makes their ISM particularly vulnerable to disruption 
by intense episodes of star formation (starbursts). Starbursts inject turbulence into the 
interstellar medium and may provide enough energy to drive an outflow of gas from the 



galaxy (Mac Low & Ferrara 1999 de Avillez & Breitschwerdt 2005 Tenorio-Tagle et al. 



2006 Cooper et al. 2008). While there have been several studies on the effects of starbursts 



on the warm (Martin 1998), hot (Ott et al. 2005b), and cold gas (Schwartz fc Martin 2004) 



in dwarf irregular galaxies, little attention has been paid to the role of magnetic fields in 
starbursting dwarf irregular galaxies. 

The presence of magnetic fields can influence the behavior of the ISM in these galaxies in 
several ways. First, the field is an important source of pressure for the interstellar medium. 



Based on observations of the ISM of the Milky Way (for a comprehensive review see Ferriere 



2001), one often assumes that the interstellar medium of galaxies is in equipartition with 



approximately equal energy in the magnetic field, cosmic rays, and turbulence. Second, 
given the shallow gravitational potential well of these galaxies, the magnetic field may be 
strong enough to be dynamically important and aid in expelling gas from the disk of the 
galaxy. The configuration of the field, however, may help or hinder the outflow. Although 



earlier work by Tomisaka (1990) and Mineshige et al. (1993) came to the conclusion that 



magnetic fields parallel to the disk could prevent the breakout of superbubbles from the 



disk, more recent 3D simulations by de Avillez & Breitschwerdt (2005) show that while a 



disk magnetic field may delay the break out of superbubbles, the amount of hot gas vented 
from superbubbles is not greatly affected by the addition of a magnetic field. For the case 



of a magnetic field perpendicular to the galactic disk, Breitschwerdt et al. (1991> 1993) 



present models of a cosmic-ray driven wind, where the magnetic field transfers momentum 



and energy from cosmic rays to thermal gas thus helping to drive the wind. Everett et al. 



(2008, 2009) later extended these models to explain the distribution of soft Galactic X-rays 



and 408 MHz emission towards the Galactic Center. The cosmic rays, acting through the 
magnetic field, may increase mass loss rates or increase the terminal velocity of the wind 



depending on the strength and geometry of the magnetic field (Everett et al. 2008, 2009). 



There are a handful of starburst galaxies with detailed magnetic field observations in- 
cluding M82 ( |Klein et al.||l~9~8~8~j |Reuter et al.||1994j ), NGC 253 ( |Carilli et al^[l992| [BeckeFaL 
1994] |Heesen et al.|[2009j ), NGC 4666 ( |Dahlem et aT]|l~997l ) , NGC 1808 ( |Dahlem et aT]|l~990| , 
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and NGC 4569 (Chyzy et al. 2006). The common features in all these objects are extended 



radio continuum emission well beyond the disk of the galaxy and spurs of polarized radio 
continuum emission extending from the disk of the galaxy to the halo of the galaxy with 
their derived magnetic field vectors roughly perpendicular to the disk. However, all these 
galaxies are quite massive - NGC 253, NGC 4666, NGC 1808, and NGC 4569 have dynamical 
masses on the order of 10 11 M , while M82 has a dynamical mass on the order of 1O 1O M . A 
low- mass galaxy like a dwarf irregular, which has a dynamical mass on the order of 10 9 M , 
provides a more rigorous test of the role of the magnetic field in a starburst. 

One of the most extreme starbursts in the local universe is the dwarf irregular galaxy 
NGC 1569. This galaxy, which is only 3.36 Mpc away (Grocholski et al. 2008), underwent 
an intense burst of star formation (star formation rate ~ 1 — 3 M yr _1 kpc~ 2 ) in the recent 
past (~ 10-40 Myr ago) QVallenari k Bomans||T996l |Greggio et al]|1998[ |Aloisi et al."l|2001 



Angeretti et al. 2005[ ). Today, NGC 1569 is in a post- starburst phase, although it is still 



vigorously forming stars (Greggio et al. 1998). The subsequent evolution and demise of the 



high-mass stars formed in the starburst has injected an enormous amount of energy into the 
interstellar medium of NGC 1569 causing massive disruption of the gas. The Ha structure 
of this highly inclined galaxy (i = 63; Stil & Israel 2002) consists of filaments of emission 



extending well into the halo of the galaxy. The most prominent Ha filament is referred to 
as the "western arm" (note that this is not a spiral arm, but the limb of an Ha bubble). 



Investigations of the Ha kinematics (Tomita et al. 1994; Heckman et al. 1995; Martin 1998 



Westmoquette et al. 2007a|b, 2008) reveal that the central regions of the galaxy are domi- 



nated by turbulent motions, while the outer regions show clear signs of a large-scale outflow 
of ionized gas. The neutral hydrogen velocity field is completely disrupted in the central 



regions of NGC 1569, but rotates as a solid-body in the outer regions (Miihle et al. 2005). 



X-ray observations ( Heckman et al.||1995 della Ceca et al.||1996 Martin et al.|2002 ; Ott et al. 
2005a[b ) have revealed the widespread presence of hot, partially shock- heated, gas bounded 
by the Ha filaments throughout the galaxy. This gas has a high a-to-Fe ratio indicating 



that the hot gas is enriched with the products of type II supernovae (Martin et al. 2002). 



Although NGC 1569 clearly has gas being expelled from its central regions and the gas likely 



has enough energy to escape (Martin 1998 Ott et al. 2005b), how much material, if any. 



will actually escape the galactic potential remains controversial. A detailed study of the 
magnetic field strength and structure of NGC 1569 is a crucial component for understanding 
the outflow of gas from this galaxy. 

The magnetic field of a galaxy can be traced using its radio continuum emission, which 
is a combination of synchrotron radiation generated by relativistic electrons spiraling along 
magnetic field lines and free-free emission generated by thermal electrons. The resolved radio 



continuum emission from NGC 1569 has been observed by Seaquist & Bignell (1976), Condon 
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( p83| , |Klem & Grave| ( p86| (1.2 cm [24.5 GHz] data only), [Israel fc de Br"uyn] ( fl988| ), and 
Lisenfeld et al. (2004). Seaquist fc Bignell| (1976) found that most of the regions they 
detected in NGC 1569 are thermal in nature; however, these observations were only sensitive 
to a limited range of spatial scales with little sensitivity to large-scale emission. |Condon 



(1983) presented VLA observations of NGC 1569 showing an extended radio continuum halo 



at 20 cm (1.4 GHz), but did not discuss them extensively. The Effelsberg 100-m observations 



of |Klein fc Grave] ( fl986| ) only resolved NGC 1569 at 1.2 cm (24.5 GHz). At this wavelength, 
they detected faint emission associated with the western arm for the first time at a frequency 



higher than 20 cm (1.4 GHz). Israel & de Bruyn (1988) presented observations of the radio 



continuum emission from NGC 1569 at 50 cm, 20 cm, and 6 cm (0.6 GHz, 1.4 GHz, and 
5.0 GHz, respectively) and discuss in detail the faint extended radio halo around this galaxy 



first found by Condon ( 1983 ) and visible in their observations at both 50 cm and 20 cm. They 



identified a break in the non-thermal radio continuum spectrum of this galaxy at 3.75 cm 



(8 ± 1 GHz). Lisenfeld et al. (2004) modeled the integrated radio spectrum (including the 
break) and concluded that the radio continuum spectrum is likely the result of cosmic-ray 
escape in a convective wind, which fits with the picture of NGC 1569 as a galaxy undergoing 
a tremendous outflow of material. 

In this paper, we present multi-wavelength radio continuum polarization observations 
of the post-starburst dwarf irregular galaxy NGC 1569, which will allow us to determine the 
strength of the magnetic field in this galaxy and trace its structure. This data set represents 
the most sensitive set of radio continuum observations of NGC 1569 to date and the first 
observations of its polarized emission at 20 cm, 13 cm, 6 cm, and 3 cm. We describe our data 
set in § [2] and present the basic properties of the total intensity and polarized continuum 
emission of NGC 1569 in § [3] In § |4| we investigate the properties of the magnetic field of 
NGC 1569. We place our observations in the context of what is known about the interstellar 
medium of NGC 1569 in § [5j Finally, in § |6j we present a summary of our work and our 
conclusions. 



2. Observations and Data Reduction 

We observed the radio continuum emission from NGC 1569 at 20, 13, 6, and 3 cm. 
The 20 and 13 cm observations were taken with the Westerbork Synthesis Radio Telescope 
(WSRT0 The WSRT is an east-west interferometer consisting of fourteen 25-m antennas. 



1 The Westerbork Synthesis Radio Telescope is operated by ASTRON (Netherlands Foundation for Re- 
search in Astronomy) with support from the Netherlands Foundation for Scientific Research (NWO). 
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Ten of the antennas are at fixed positions, while the other four antennas may be moved along 
an east-west track. At 6 and 3 cm, we used both the Very Large Array (VLAj^J which is an 
interferometer with twenty-seven 25-m antennas arranged in a Y-shape, and the Effelsberg 
100-m Radio Telescopd^ which is a single-dish telescope. The WSRT and VLA observations 
provide us with the spatial resolution necessary to detect small-scale magnetic field structure 
that would be depolarized by the large beam of the Effelsberg observations. The Effelsberg 
observations are necessary to recover the emission on the largest spatial scales, which may 
be resolved out by the VLA at 6 and 3 cm. At 20 and 13 cm, the largest angular scales 
imaged by the WSRT are approximately 19' and 17', which are both much greater than the 
size of the galaxy, so there is very little chance of large-scale emission from NGC 1569 being 
resolved out at these wavelengths. Below we discuss the calibration and imaging of the data 
in detail. 



2.1. 20 cm Data 

The 20 cm data were taken with the WSRT as part of project R01A/12 (PI: S. Miihle) on 
5-6 and 6-7 November 2000. The project also included observations of the radio continuum 
of NGC 1569 at 49 cm, but the 49 cm data were unusable due to extensive radio frequency 
interference. To calibrate the instrumental polarization, each observing session included an 
observation of a source with a known polarization angle and an unpolarized source. These 
sources were also used to calibrate the flux density scale. The time on-source for each day 
was 12 hours. To maximize the coverage of the u-v plane, we used two different array 
configurations: the 72 m configuration and the 36 m configuration. The DCB backend was 
configured to produce eight lOMHz-wide intermediate frequency (IF) channels each with full 
polarization products. 

Table [T] summarizes the 20 cm observations. The data were calibrated and imaged using 



the NEWSTAR package (Brouw 1996 ), following standard calibration procedures and using a 
hybrid cleaning method, which includes self-calibration of the gains and phases of the Stokes 
I, U, and V data. In Stokes Q, only the phases were self-calibrated, since a self- calibration of 
the gains would force Stokes Q equal to zero in the case of linear feeds (Q=XX-YY). Table [2] 
gives details for the final images. The upper left panel of Figure [T] shows the total intensity 



2 The VLA is part of the National Radio Astronomy Observatory. The National Radio Astronomy Obser- 
vatory is a facility of the National Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 

3 Operated by MPIfR (Max-Planck-Institut fur Radioastronomie) . 



-7- 



map of the 20 cm data. The final image includes low-level (20 /xJy beam 1 ) ripples in the 
total intensity that are artifacts of the imaging process for this data. 

2.2. 13 cm Data 

The 13 cm data were obtained with the WSRT (project ID: R03A/25, PI: S. Miihle) 
on 4-5 December 2002 and 18-19 January 2003. The time on-source for each observing 
session was 12 hours. Observations of a polarized source and an unpolarized source were 
taken before and after each session to calibrate the instrumental polarization as well as the 
bandpasses and the flux density scale. A summary of the WSRT 13 cm observations used in 
this paper and their calibration is given in Table [3j The data were taken in the 36-m and 
72-m configurations and DZB backend was configured to produce eight 20 MHz-wide IFs, 
each with full polarization products and 64 channels. 

Since the 13 cm receivers on the WSRT have circular feeds, we calibrated the data using 
the standard calibration routines in AIPS, which, in general, assume circularly polarized 
feeds. Reducing WSRT data in AIPS has two important differences compared to reducing 
typical VLA observations. First, the WSRT is a redundantly spaced array, so observing phase 
calibrators periodically during the observing session is not necessary because the redundant 
baselines provide multiple measurements of the same visibility. Self-calibration is used to 
solve for the variation of phase as a function of time. The second complication is that the 
WSRT dishes are equatorially mounted, so the parallactic angle does not change as a function 
of time. Therefore, a single observation of an unpolarized source must be used to calibrate the 
instrumental polarization. This is in contrast to the VLA where the polarization response of 
the antennas is usually calibrated using a weakly polarized source observed over a large range 
in parallactic angle. Practically, it means that the task LPCAL must be used to calibrate 
the instrumental polarization instead of the task PCAL. 

We used the AIPS task CALIB to self-calibrate the data using models of the galaxy 
generated by imaging and cleaning the data. We did several iterations of phase-only self- 
calibration with solution intervals of between 5 and 10 minutes and restricted the solutions 
to those baselines where the signal-to-noise ratio was greater than 3. After the phases had 
been improved, we did a final round of amplitude and phase self-calibration with our best 
model of the source, which included as much of the diffuse flux as possible. The solution 
interval was between 20 and 30 minutes and only baselines where the signal-to-noise ratio 
was greater than 5 were included in the solution. Each observation period was self-calibrated 
separately and the individual periods combined together. Finally, the combined data was 
amplitude and phase self-calibrated with a solution interval of 12 hours including only those 
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baselines with signal-to-noise ratio greater than 5 to ensure a uniform amplitude scale. 

We produced two final images: a robust=0 image and an image matched to the resolution 
of the high-resolution 20 cm image (~ 13"). The robust=0 image has a weighting between 
uniform and natural and provides the best compromise between the lower noise given by 
natural weighting and the minimization of the sidelobes by uniform weighting. Table [2] lists 
the details of the final images. The upper right panel of Figure [T] shows the robust=0 image. 



2.3. 6 cm & 3 cm Data 

The interferometer data at these wavelengths were taken with the VLA under program 
IDs AM643 (D array observations) and AM694 (C array observations). The PI for both 
programs was S. Miihle. A full set of polarization products (RR, RL, LR, LL) were produced 
for each of 2 IFs with a bandwidth per IF of 50 MHz. The total time on-source was 3.97 
hours for the 6 cm (4.8 GHz) D-array observations, 14.7 hours for the 6 cm (4.8 GHz) C-array 
observations, and 10.5 hours for the 3 cm (8.4 GHz) D-array observations. We lost 2 hours 
of time on source during the 3 cm observations on 23 September 2000 due to high winds. 
We calibrated the data in AIPS following the standard calibration procedures detailed in 



the AIPS cookbook (Greisen 2006). We self-calibrated the data in AIPS using the same 



procedure as for the 13 cm data (see § 2.2). Table [4] summarizes our observations and 
calibration parameters. 

Whenever possible, we used 1331+305 (3C286) as our absolute polarization angle cali- 
brator, since the polarization angles of all other calibrators are derived from this source. If 
1331+305 was not available, 0521+166 was used as an absolute polarization angle calibrator. 
The polarization angle of 0521+166 in the last row of Table|4]was determined using observing 
sessions where we observed both 0521+166 and 1331+305. The absolute polarization angle 
for 0521+166 on 9 August was determined from observations on 10 August, on 11 and 17 
August from observations on 19 September, on 23 September from observations on that day. 
Our observations on 19 September indicate that the absolute polarization angle of 0521+166 
varied by ±2° over the entire observing run. 

To correct the 6 cm and 3 cm VLA data for the effects of missing short spacing^} we 
obtained 6 cm (4.85 GHz) and 3 cm (10.45 GHz) data with the Effelsberg 100-m radio tele- 
scope in December 2000 and January 2001, respectively. The data at both wavelengths were 



4 It is not necessary to correct the 20 cm and 13 cm data for the effect of missing short spacings since the 
largest angular scale imaged by the WSRT at these frequencies is much larger than the size of the galaxy. 



-9- 



reduced and calibrated using the N0D2 data reduction package (Haslam 1974). The 6 cm 



(4.85 GHz) Effelsberg receiver has two feeds, which allows one to eliminate atmospheric noise 



using the dual-beam technique described in Emerson et al. (1979). The telescope scans in 



azimuth to align the two feeds. To remove raster noise, we averaged together individual 
maps taken across a range of parallactic angles. Atmospheric noise at 3 cm (10.45 GHz) 
was minimized by using the software beam-switching technique on the data from the four 
feeds of this receiver. The residual scanning noise in the individual images was reduced by 



down- weighting the regions in Fourier space where raster noise is present (Emerson & Grave 



1988). The individual maps were then averaged to produce the final maps. Table 5 gives a 



synopsis of our Effelsberg observations. 

Since the VLA and Effelsberg 3 cm images were taken at two different frequencies (8.46 
GHz and 10.45 GHz), we scaled the Effelsberg Stokes I image by the spectral index deter- 
mined by Lisenfeld et al. (2004): a = —0.47. We did not scale the Effelsberg Stokes Q and 
U images because we do not know how NGC 1569's polarized emission changes between 8.46 
GHz and 10.45 GHz. We can, however, constrain the rotation measure between these two 
frequencies. For a change of 10° in the angle of the polarization vector, there would have to 



be a rotation measure of at least 400 rad m 2 . In 5 4.2, we determine the rotation measures 



for NGC 1569 between 6 cm and 3 cm. Most of the measured values are less than 400 rad m 2 



(see Figure 10), so not scaling the flux in the Effelsberg 10.45 GHz Stokes Q and U images 
should not be problematic. However, rotation measures greater than, or comparable, to 400 
rad m~ 2 may be affected by combining unsealed Stokes Q and U data taken at two different 
frequencies. 

To combine the interferometric and single-dish data, we used the immerge task in Miriad, 
which combines the images in the spatial frequency domain. In this approach, the different 
calibration scales of the Effelsberg and VLA data are dealt with by scaling the low-resolution 
visibilities so that they match the high-resolution visibilities in the region of the u — v plane 
where the observations overlap (0.4 to 1.6 kA for the 6 cm data and 0.725 to 2.8 kA for 
the 3cm data). The scaling factors for the visibilities in the overlap region as determined 
by immerge were 0.8196 (6cm, robust=0 image), 0.856 (6cm, 13" image), 0.9263 (3cm, 
robust=0 image), and 0.9330 (3cm, 13" image). We used the scaling factors determined by 
the Stokes I images for the Stokes Q and U images. 

Details on the final 6 cm and 3 cm images are given in Table [2j Figures [T] and [2] show 
the robust=0 weighted 6 cm and 3 cm images. The total fluxes of the final images are within 
10% of the total fluxes measured with Effelsberg (see Figure [3]). 



3. Properties of the Radio Continuum Emission 



The total intensity of the radio continuum emission from a source gives one information 
about the combined free-free and synchrotron emission in that source. By separating the 
synchrotron emission from the free-free emission, one can estimate the total magnetic field 
strength (see §£4.1 and 4.2). Since free-free emission is unpolarized and synchrotron emission 
is polarized, any polarized emission is inherently synchrotron emission and traces of the 
uniform component of the magnetic field (see § 4.3). 



3.1. Total Intensity 



Previous resolved studies of NGC 1569 at radio continuum wavelengths include Seaquist 



fc Bignell (1976), Condon (1983), Klein & Grave (1986) (1.2cm [24.5 GHz] data only), Is 



rael & de Bruyn (1988), and Lisenfeld et al. (2004). The details of these observations are 



shown in Table [6] where necessary, we converted sensitivities in brightness temperatures to 
sensitivities in mJy beam -1 using the given frequencies and beam parameters. Our obser- 
vations reach surface brightnesses two (20 cm observations) to eighty (13 cm observations) 
times fainter than previous observations. Although our resolution is not as high as some pre- 
vious observations, we do have better u-v plane coverage, making our observations sensitive 
to both compact sources and extended emission. 

We have detected radio continuum emission from NGC 1569 at all four observed wave- 
lengths. In Figure |3j we plot the total flux detected by our measurements along with radio 
continuum flux measurements collected from the literature by Israel & de Bruyn (1988) and 



Lisenfeld et al. (2004). Our measured fluxes at 20cm (1.4 GHz) and 13cm (2.7 GHz) agree 
well with the measured fluxes from the literature, justifying our assumption not to correct 
the WSRT observations at these frequencies for the effects of missing short-spacings. The 
6 cm (4.86 GHz) and 3 cm (8.46 GHz) fluxes from our combined VLA and Effelsberg images 
are -9.5% (6 cm) and -8.8% (3 cm) less than the measured Effelsberg fluxes. (Note that 
the low point from the literature at 3 cm was an interferometer measurement that did not 
correct for missing short-spacings.) All our observations are at wavelengths greater than or 
approximately equal to the break in the radio continuum spectrum at approximately 3.75 cm 



8 GHz) first noted by Israel & de Bruyn (1988) and modeled by Lisenfeld et al. (2004). 



Our total intensity images at all four frequencies are shown in Figures [T] and [2] Note 
that we show the highest resolution image at each frequency, so the beam size varies between 
the images. At 20 cm, we detect both the bright central regions and the extended diffuse 



component first seen by Condon (1983) and first discussed by Israel & de Bruyn (1988). 
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There are two bright peaks of emission in the central regions of the galaxy. The brightest 
peak is on the northwestern side of the disk and the fainter, slightly elongated peak is on the 
southeastern side of the disk. In our 6 cm, robust=0 image, which has the highest resolution 
of all our images, we resolve the southeastern peak into two sources. Our observations 



increase the extent of the diffuse component of the galaxy from ~ 3' (Israel & de Bruyn 



1988) to ~ 4'. This corresponds to a physical extent of approximately 4.4 kpc off the 



plane of the disk (assuming an inclination of 63°; Stil & Israel 2002), which is extraordinary 
considering the disk is only about 3' (2.9 kpc) long in the radio continuum. We also detect 
the western arm of the galaxy at all observed frequencies, where previous observations only 
saw it at 20 cmj^] As one goes from 20 cm to 3 cm, the morphology of the diffuse component 
of the radio continuum emission changes from box-shaped to more oval-shaped. At shorter 
wavelengths, the western arm becomes increasingly prominent and the extent of the radio 
continuum halo decreases. 

Spectral index maps of the observed total intensity radio continuum emission can provide 
an important diagnostic of the processes producing the radio continuum emission. We define 
the spectral index, a, as S oc u a , where S is the source flux and v is the frequency; thermal 
emission has a spectral index of —0.1 and synchrotron emission has a spectral index of 
~ —0.7. By looking at the observed spectral index, we can get an idea of whether thermal 
or synchrotron emission dominates in a particular region. 

Spectral index maps, however, are significantly affected by small systematic changes in 
flux. For example, a 5% change in the flux level in each total intensity map going into the 
spectral image can change the spectral index by as much as 0.18. Because of the artifacts 
remaining in the 20 cm and 13 cm images and the systematic uncertainties in creating spectral 
index maps from data from two different telescopes (especially ones as different as the VLA 
and the WSRT), we created a spectral index map for NGC 1569 based only on the 6cm 
and 3 cm data. To make sure that both images sampled the same spatial scale, we used 
only the interferometric data and restricted the u-v range of the data imaged to 0.74 to 21 
kA; this range is the u-v range common to the observations at 13 cm, 6 cm, and 3cmJ^] The 
beam sizes in each image were slightly different because differences in the sampling of the 
u-v plane in each image. We smoothed both images to the resolution of the image with the 
largest beam (9.34" by 7.54"). 

Figure [4] shows the spectral indices between the 6 cm and 3 cm data. Only spectral 



5 Klein & Grave ( 1986 ) saw a hint of this feature in 24.5 GHz single-dish data. 

6 The range common to the 6 cm and 3 cm is 0.74 to 29.3 kA, which would increase the resolution slightly. 
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indices where the formal error on a was less than 0.2 were used. The formal error on a is 



log(z/i/V 2 )_ 



Si 



02 

s 2 



where S, u, and a are the flux in a pixel, the frequency, and the noise for a particular 
image and the subscripts denote the first and second images. The values in our final spectral 
index map agree with those in a spectral index map created using data that was not self- 
calibrated, showing that our self-calibration of the 6 cm and 3 cm data has not introduced 
major uncertainties into our spectral index maps. 

The spectral indices in the galactic disk are close to thermal with flattest spectral 
indices offset by 14" (a little more than a beam) to the north from the main ridge of 3 cm 
continuum emission. Since the southern half of the galaxy is pointed toward the observer 



and the northern half pointed away (Martin et al. 2002), this spectral index asymmetry 



about the major axis of the disk can be explained as free-free absorption of synchrotron 
emission from behind the disk by the disk itself. Free-free absorption flattens spectral indices 



(Condon 1992). In general, as one moves along the minor axis of the galaxy, away from the 
disk, the spectral indices become steeper, eventually reaching values of approximately -1 
at the southern edge of the galaxy. There are, however, two filamentary regions of flat 
spectral indices extending to the outer edge of the spectral index map. These flat spectrum 
filaments are aligned with the outer edges of an Ha bubble (Bubble B in the terminology of 
Westmoquette et al.||2008 , see also Figure [7]). The spectral index trends seen here agree with 



the general trends seen in the higher resolution data of Lisenfeld et al. (2004) supporting 



their theory that the break in the total non-thermal emission spectrum of NGC 1569 is due 
to the convective transport of cosmic rays into the halo. 



3.2. Polarized Emission 

The observations presented in Figures [5] and [6] are the first published observations of the 
polarized emission of NGC 1569. We detect polarization at 20 cm (36 to 80 /zJy beam -1 ), 
13cm (49 /iJy beam -1 to 83 /iJy beam -1 ), 6 cm (30 to 108 /iJy beam -1 ), and 3cm (31 
to 137 /iJy beam -1 ). At 3 cm and 6 cm, there is polarized emission over a large portion 
of the galaxy, while at 20 cm the polarized emission is much spottier and at 13 cm there 
are only a few small patchies of polarized emission. Note that the strongest polarization 
detections at 6 and 3cm were with the images with 13" resolution, i.e., similar resolutions 
to the 20 cm data. Figures [5] and M show the distribution and direction (E- vectors) of the 



13 



polarized intensity at the observed wavelengths^ To correct for the bias in the estimated 
polarized intensity, we use the POLCO task in AIPS with a signal-to-noise cut of three. This 
signal-to-noise cut ensures that all polarization vectors are significant at the 99% confidence 



level (Simmons & Stewart 1985; Vaillancourt 2006). Decreasing the signal-to-noise cut to 



two did not add appreciably more vectors. Despite NGC 1569's position near the Galactic 
plane, (I, b) = (143?7, 11?2), there is fortunately very little foreground Faraday rotation along 



this line of sight. From the Galactic Faraday rotation map of Johnston- Hollitt et al. (2004), 
we estimate that the foreground Faraday rotation is approximately +10 rad m -2 . Examining 
the polarized sources in our 20 cm and 13 cm data (unfortunately, the only polarized source 
in our 6cm and 3 cm data is NGC 1569 itself), we find an average rotation measure between 
the 20 cm and 13 cm data of -8.7 rad m -2 with a standard deviation of 36.7 rad m~ 2 , which 



agrees with the value from the Johnston- Hollitt et al. (2004) map within the errors 



The polarized emission at 20 cm is shown in the left hand panel of Figure |5j The 
percent polarization in this image ranges from 2% at the center to 70% at the edges. The 
central region of NGC 1569 is bright enough to have significant polarized emission resulting 
from residual instrumental leakages rather than emission intrinsic to the galaxy. Based on 
measurements of unpolarized sources in the field, we estimate this leakage at approximately 
1%. Therefore, we require the polarized emission at 20 cm to have a fractional polarization 
of greater than 1% as well as a signal-to-noise ratio of three to be considered significant. The 
emission in the final image could be depolarized if the polarization vectors rotate significantly 
over the observing frequencies that are combined to produce the final image; this effect is 
known as bandwidth depolarization. To determine if the 20 cm data is affected, we imaged 
the first and last intermediate frequencies (IFs) separately (see Table [I] for a summary of our 
observing setup). These two IFs are each 10 MHz wide and 70 MHz apart. Comparing the 
polarization vectors for the two IFs, we do not see any rotation of the polarization vectors 
greater than a few degrees. Therefore, creating the final image by combining the data from 
all the IFs has not resulted in significant depolarization. 

The major polarization features at 20 cm are associated with Ha features in the galaxy. 
Figure [7] shows the positions of the various Ha bubbles that have been identified in the 
literature. The Ha emission is included in both Figure [5] and [7] as a point of comparison 
between the figures. There are polarization features on the northwest and southwest edges 
of the disk which correspond to the bases of bubbles B and A. Moving northwest away from 
the disk there are regions of significant polarized emission that outline the edges of bubble 
B as well as some that may be associated with the region of Ha filaments just to the east of 



7 We plot E- vectors in Figures[s]and[6]because they are the quantity actually measured by the observations. 



We rotate the E-vectors by 90° to estimate the direction of the magnetic field in Figure 10 
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this feature. We also see significant polarized emission outlining the bubble A. 

The 13 cm emission is shown in the right panel in Figure [5] The polarization ranges 
from 1.4% for vectors near the center of the image to 42.1% for the clump of significant 
vectors in the northwestern portion of the galaxy. We used one of the calibration sources 
(3C147) to estimate the residual instrumental leakages at 13 cm as 0.5% and eliminated 
all vectors with smaller percent polarizations. This cut affected only one vector located 
on the brightest peak of 13 cm emission. There is very little significant polarization in the 
13 cm data besides the polarization vectors in the northwestern portion of NGC 1569 that 
correspond to locations of significant polarization vectors at 20 cm. The lack of significant 
polarized emission is not surprising since the 13 cm has the highest noise level of all our 
data. When the 20 cm polarized intensity is calculated assuming a noise cutoff similar to 
that of the 13 cm data we see little polarization, which confirms that the lack of polarization 
at 13 cm is most likely the result of the relatively poor sensitivity of this data. To make sure 
this is not the result of bandwidth depolarization, we have imaged first and last 20 MHz 
IFs by themselves and compared polarized emission in the two images. Neither one of the 
single IF images shows any significant polarization, so significant polarization in one IFs is 
not canceling out significant polarization in another IF. 



The most polarized emission is seen in the 6 cm and 3 cm data (Figure [6]). The percent 
polarization seen in the 6 cm and 3 cm data is well above the estimated leakage of the VLA 



<1% near the center of the field; Ulvestad & Taylor 2007). Therefore, the only cut we have 



applied to this data is to require the signal-to-noise ratio for the polarized emission to be 
greater than three. At both 6 cm and 3 cm, there is no visible rotation of the polarization 
vectors between the two IFs, so bandwidth depolarization does not seem to be an issue. 
There are some significant polarization vectors outside the 3a contour line for the 3 cm total 
intensity. For the 3 cm data with the 13" resolution, however, the error in the Stokes I image 
is larger than estimates of the theoretical noise in Stokes I based on the noise in the Stokes 
Q and U images. Therefore, we are relatively more sensitive to polarized emission than to 
total intensity emission. 

In general, the polarized emission at 3 cm is strong and present over a large fraction of 
the galaxy. The polarized emission at this wavelength goes from 1-2% at the center to up to 
40% and greater on the western edge and 10-30% in the halo of the galaxy. There is little 
and/or weak polarized emission coincident with the western Ha arm and another break in 
the emission at approximately 04 h 30 m 45 s , 64°50'20". The central regions of the 3 cm data 
appear to be largely unpolarized except for some polarized emission associated with the 
main radio continuum emission peak. This depolarization could be due either to Faraday 
depolarization or to an inherently turbulent field. We discuss these possibilities further in 
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§ [5] in the context of NGC 1569's Ha kinematics. The majority of the polarized emission 
is in the southern half of the galaxy, which is unsurprising given that X-ray observations 
have found that the southern portion of the galaxy is inclined towards us and the northern 



part is inclined away (see Figure 13 in Martin et al. 2002). The polarized emission from the 



northern half of the galaxy has a longer path length through ionized gas and thus is subject 
to more internal depolarization. There is particularly strong polarization on the western 
edge of the Ha arm. This feature is also seen in higher resolution images not included here. 
The region on the interior of the Ha arm (Bubble A in Figure [7]) is polarized as well as the 
base of bubble F to the southeast of the disk. In the northern half of the galaxy, there is 
polarized emission on the northeastern edge corresponding to Bubble E and to the north of 
the main radio continuum peak, which corresponds to Bubble B. 

The polarized emission at 6 cm is patchier than at 3 cm and ranges from 1-2% in the 
inner regions of the galaxy to 15-20% in the outer regions of the galaxy to greater than 40% 
in the polarized arm. Interestingly, the E- vectors in common between the 3 cm and 6 cm data 
show little sign of Faraday rotation except for the polarization in the northern half of the 



galaxy where there is significant rotation (see § 4.3). There is significantly more polarized 
emission north of the plane of the galaxy at 6 cm than at 3 cm. Since emission at longer 
wavelengths tends to trace emission at smaller Faraday depths, this emission is probably 
closer to the observer along the line of sight than the emission at 3 cm. At 6 cm there, is 
still polarized emission surrounding the Ha arm as well as polarized emission at the base of 
Bubble F in a region of extensive Ha filaments and in the Bubble E region. These regions 
of emission, however, are less extended than those seen in the 3 cm data. The central region 
of NGC 1569 at 6 cm has no significant polarization. 

In short, our observations reveal that more polarized emission is seen in NGC 1569 at 
shorter wavelengths than at longer wavelengths and that the polarized emission is associated 
with Ha bubbles and filaments in the galaxy. The first result is common to radio continuum 
polarization observations of other galaxies ( Beck et al.||1996 ). Depolarization due to Faraday 



dispersion and differential Faraday rotation increases with increasing wavelength, so radio 
continuum polarization observations of galaxies at longer wavelengths show less polarization. 



We defer discussion of the second result to 5 4.3 
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4. Magnetic Field Structure and Strength 

4.1. Estimating the Synchrotron Flux 

To determine the magnetic field strength in NGC 1569 and thus assess the relative 
importance of the field to the ISM, we must estimate the portion of the total radio continuum 
flux due to synchrotron emission and the portion due to free-free emission. To do this, we 
estimated the thermal component of the radio continuum emission using the Ha emission 
from the galaxy. We used a flux calibrated Ha image kindly provided by Deidre Hunter 



(Hunter & Elmegreen 2004) to determine the Ha flux in each pixel. We corrected the 



resulting Ha flux for Galactic extinction using a value of An a = 1.26 (Relano et al. 2006). 



Since the internal extinction of NGC 1569 is small (A^a < 0.6 mag) and largely confined 



to the central regions of the galaxy (Relano et al. 2006), we do not attempt to correct for 



internal extinction. Not correcting for internal extinction will cause us to underestimate the 
Ha flux by, at most, a factor of 1.74 in the regions with the greatest extinction, which are 
the region surrounding the central peak and the region connecting the two brightest radio 



continuum peaks (see Figure 9 of Relano et al. 2006). Note that we do not need to take 



into account non-photoionization processes since the calculation of the conversion between 
Ha emission and free-free emission does not depend on the ionization source. 

We calculated the expected thermal emission at our observing frequencies from the given 
Ha flux using the expressions given in Appendix A of Hunt et al. ( 2004) j^] which explicitly 
include the contribution of ionized helium: 



fu 



1.16 1 + 



n(#e H 



T 



10 4 K 



0.617 



GHz 



-o.i 



Ha 



10 -12 erg cm -2 s 1 



(2) 



mJy J \ n(H + 

In Equation (j2j), f v is the radio continuum flux, n(He + )/n(H + ) is the ratio of ionized 
helium to ionized hydrogen, T is the temperature of the region, v is the frequency of the 
radio continuum observations, and fua is the Galactic-extinction-corrected Ha flux. We use 



a n(He + )/n(H + ) ratio of 0.087 (Martin & Kennicutt 1997). The expressions from Hunt 



et al. (2004| are valid for temperatures between 10,000K and 20,000K and densities between 
100 cm -3 and 1000 cm -3 . We assume a temperature of 10,000K, which is typical of HII 



regions. This temperature is 10% lower than the lowest temperature that Devost et al. (1997) 
was able to measure for selected regions of NGC 1569. If we increase the temperature to 
11,100 K, this increases the estimated thermal emission by only 7%, since the temperature 
dependence in this equation goes as T 0617 . Westmoquette et al. (2008) shows that the 



densities of ionized gas in NGC 1569 range from 100 cm (the lowest density probed by the 



See also 



Caplan & Deharveng 



(1986 


) and 


Condon 


(1992 



-17- 



their diagnostic) to 1000 cm 3 . We smoothed the result of Equation §2§ to the resolution of 
the observations using the AIPS task CONVL. The synchrotron flux is 

Fnth = F to t [ 1 — ] (3) 

V ftot J 

where F nt h is the flux of the synchrotron (non-thermal) emission, F tot is the total flux ob- 
served, and F t h is the thermal flux calculated using Equation §2§. 

We plot the estimated thermal fraction (F th /F tot ) in Figure [8] Note that we set thermal 
fractions greater than one to one and thermal fractions less than zero to zero. This does 
not significantly change the resulting images and only affects regions where the emission is 
faint, i.e., the edges of the map. As one would expect, the total emission at both 20cm 
and 13cm is dominated by synchrotron emission, i.e., has a low thermal fraction, with the 
maximum thermal fraction (40%) southeast of the location of the radio continuum peak. 
The halo of the galaxy has significantly more synchrotron (less thermal) emission than the 
disk. As we go to shorter wavelengths, the amount of thermal emission increases to 50% 
in the region southeast of the radio continuum peak. Note that the inner portion of the 
western arm has a significant amount of thermal emission, while the outer portion is largely 
synchrotron emission. The outline of Bubble A appears very clearly in the 3 cm map of 
the thermal emission, which provides additional confirmation that we are able to accurately 
determine the location of the thermal flux with the procedure described earlier in this section. 
The entirely thermal regions on the northern and eastern edges of the 3 cm map are due to 
the 3 cm diffuse emission being slightly underestimated there. Although we have been very 
careful about reconstructing most of the diffuse flux, it is impossible to reconstruct it all, 
which results in the thermal fraction being overestimated in these regions. 



4.2. Magnetic Field Strength 

The intensity of the non-thermal (synchrotron) emission can be used to estimate the 
magnetic field strength in a galaxy. The intensity of synchrotron emission at a particular 
frequency goes as 

I v oc lN Q B x ~ a v a (4) 

where / is the line of sight through the source, N is the scaling factor for the electron energy 
spectrum (N(E) dE = NoE~ 7 dE, 7 = 1 — 2a), B± is the magnetic field strength in the plane 
of the sky, and v is the observing frequency (Ginzburg & Syrovatskii 1965J ) p] For a par- 



ticular value of I u only the quantity N B) a can be estimated, so a particular synchrotron 



3 We use a sign convention for a opposite to the one typically used. 
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intensity may correspond to a large electron density and a small magnetic field or vice versa. 
To disentangle the magnetic field strength, one traditionally employs the minimum energy 
assumption where the total energy (the energy of the magnetic field plus the energy of 
the cosmic ray electrons plus the energy of the cosmic-ray protons) is minimized. For the 
typical range of 7, the minimum energy is similar to the energy derived assuming equipar- 



tition between cosmic rays and the magnetic field (Beck & Krause 2005). The minimum 



energy/equipartition calculation requires a further assumption that the energy spectrum of 
the cosmic-ray protons is the energy spectrum of the cosmic-ray electrons scaled by a con- 
stant. However, this is not true in the Milky Way, where the electron spectrum is steeper 



than the proton spectrum (see Figure 1 in Hillas 2006). Beck & Krause (2005) present a 



revised equipartition estimate of the magnetic field strength that uses the number density 
ratio of electrons and protons rather than the more typical energy density ratio to estimate 



the field strength. They use the equation below (Equation 3 in Beck & Krause 2005) to 
estimate the total equipartition magnetic field strength: 



4tt (1 - 2a) (go + 1) U E l p +2a (v/2 Cl ) 
(—2a — 1) c 2 I cAi) 



l/(3-«) 



(5) 



where B t is the total equipartition magnetic field in G, Kq is the number density ratio of 
protons to electrons, E p is the rest energy of the proton, c\ is a constant equal to 6.26428 x 
10 18 



erg 



2 s" 1 G" 



C2 is a constant tabulated on page 232 of Pacholczyk (1970) as C5, I is 



the path length through the galaxy, and c±(i) is a function correcting for the inclination of 
the region with respect to the sky p^j [Beck &: Krause (2005) caution that this formula may 
not be strictly applicable to starburst galaxies, since synchrotron losses are important in 
high-density regions, and that detailed modeling is needed to determine the true magnetic 



field of a starburst. However, based on the work of Thompson et al. (2006), we believe that 



our magnetic field strength estimates for NGC 1569 will not be significantly underestimated. 



Thompson et al. (2006) estimate, for a large number of galaxies from normal spirals to 



ULIRGs, what they call an "equipartition" magnetic field using the assumption that the 
magnetic energy density is comparable to the gravitational pressure of the ISM and compare 
it to a minimum energy magnetic field estimated using a method similar to Equation ([5]). 
The Thompson et al. "equipartition" field is the largest magnetic field a galaxy can have 
and still be gravitationally stable object. In Thompson et al. Figure 1, NGC 1569 lies on 
the locus where both the "equipartition" magnetic field and the minimum energy magnetic 
field are similar, far from the position of galaxies like NGC 253, M82, and Arp 220, whose 
minimum energy magnetic field estimates are well below the Thompson et al. "equipartition" 
field. 



10 We use the opposite sign convention for a as 



Beck & Krause 



(2005) 
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We use the method of "magnetic maps" described in Chyzy (2008 ) to make a map of the 
total (B t ), uniform (B u ), and random (B r ) magnetic field strengths as functions of position 



in the galaxy along with the revised equipartition equations given in Beck & Krause (2005) 



to calculate the magnetic field strength. Equation ([5| assumes that the magnetic field is 
either all uniform or all random. Since the magnetic field of NGC 1569 is clearly a mixture 
of the two in many regions, we would like to derive a version of Equation (|5j that accounts 
for both the uniform and random field components. Equation ^ is derived from the more 
general Equation A. 16 in Beck & Krause (2005), which gives the cosmic ray energy density. 
By equating the magnetic field energy density (e# = B 2 /8n) to the cosmic ray energy density 
{^cr)i one obtains 



B[ 

8tt 
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(6) 
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Note that a here equals 7 = 1 — 2a in Beck & Krause (2005) and that c 2 is a function of a. 
Since Equation [8] includes B t , B±, and /, it not only needs to take into account the geometry 
of NGC 1569 and its associated magnetic field, but also the contribution to the total field 
from both the uniform and random components. We model the geometry of NGC 1569 as 
an infinite cylinder rotated on the plane of the sky and inclined with respect to the line of 



sight. Appendix A.l gives an equation for the line of sight distance through the galaxy. The 
geometrical correction to obtain the total field from the observed field is more complex. For 
an isotropic random field, the correction factor (c 4 ) is 



C4 



B r 



B r 



(l-a) 



(l-a)/2 



(9) 



For the regions of NGC 1569 where we do not see any polarization, we assume that the 
field is completely random. In this case, the field strength is given by Equation ^ with the 
geometrical correction factor given by Equation ([9]). For a completely uniform field, 



c 4 
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(10) 



where PAb is the position angle of the magnetic field on the sky. For the rest of this paper, 
C4 will refer to the version for a completely uniform field. For a derivation of this quantity, 
see Appendix A. 2 
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We cannot assume that regions of polarized emission are completely uniform; there will 
still be a contribution of the random component in these regions. To deal with this, we first 
need to estimate the ratio, q, of the random magnetic field (B r ) to the uniform magnetic 
field in the plane of the sky (B Ut ±). To estimate the ratio q, we use the observed percent 
polarization (p) and the relation from Sokoloff et al. ( |1998 ) 



V = Po 



1 + (7/% 2 



1 + 3q 2 + (10/9)g 4 



where po is the intrinsic polarization of synchrotron emission (p = (3 — 3a)/ (5 — 3o;))p] 



Equation (11) is derived assuming that the cosmic rays are coupled to the magnetic field 



(Sokoloff et al. 1998 Beck et al. 2003). The more commonly used version of this formula 



(Burn 1966) 



V = Po 



2„2 



1 + 



12) 



does not include this assumption. Both the Equations (11) and (12) assume isotropic fluc- 



tuations. Using q, one can now recast the total magnetic field B t in terms of the B r and B u . 



Using the relations B 



B 2 r + Bl B] 



B 2 ± + B\ _|_ and the geometrical parameters given 



in Equations \5\ and 10 
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(16) 
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Substituting the above expressions for Bl and B^_ ^ into Equation (8), we obtain an ex- 
pression giving uniform field strength in a region where we have both a uniform and random 
field 



B u 
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Beck & Krause 



(2005 )'s version of this equation has a sign error, cf. 



Ginzburg & Syrovatskii 
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In the regions where we have both a random and uniform component to the field, we calculate 



q using Equation (11) and then calculate B u from Equation (19) using the total synchrotron 



emission and the value of C4 determined from the polarization angle of the emission using 



Equation (10). The total and random fields can then be derived from the values of B u , q, 
C4, and a. 



B r 
B t 



q~DuC 4 



(Bt + Bi) 



(20) 
(21) 



For the other parameters in Equations 19 20 and 21 , we determined the position angle 
and center of the disk by fitting ellipses to the total intensity radio continuum image at 



3 cm using the IRAF (Tody 1986, 1993) package ellipse. The position angle of the disk is 



-62.1°, which corresponds to an angle (/3) between the plane of the disk and the RA axis of 
27.9°, and the center of the galaxy is (4 h 30 m 49!25, 64°50'53'/0). The diameter of the disk of 
NGC 1569 was estimated from the 3 cm image to be 2.1 kpc. We assume that the distance 



and inclination of the galaxy are 3.36 Mpc (Grocholski et al. 2008) and 63° (Stil & Israel 



2002), respectively. Since even small errors in the determination of the non-thermal fluxes 



of individual pixels will correspond to large errors in aj^] we estimate a value for a by 
calculating a for the total non-thermal flux at 6 cm and 3 cm. This method gives a value of 
a = —0.97. We take K to be 100, which is a typical value in the solar neighborhood, but 



which may be as low as 40 in regions of strong shocks (Beck & Krause 2005). Combining 



all these values, we finally determine B u , B r , and B t for each pixel in the 3 cm maps using 



Equations 19, 20, and 21 



Our magnetic maps based on the 3 cm images with 13" resolution are shown in Figure [9j 
The maximum total magnetic field strength corresponds to the peak in radio continuum 
emission at 3 cm and reaches a value of 38 /iG. In the extended halo of the galaxy, the total 
magnetic field has a value of 10-15 /iG. Examining the plot of q values, i.e., the ratio of 
the random field to the uniform field in the plane of the sky, we see that q is high (~ 20) 
in the central regions of the galaxy, suggesting that the random field dominates in these 
regions, although there could be significant Faraday depolarization in this region due to the 
high electron density. In the outer edges, however, q is only 1-2, so the regular field is quite 



strong in these regions. This effect is also seen in the magnetic fields of spiral galaxies (Heiles 



1996). The regular field in the synchrotron arm ranges between 9 /zG at the southern tip to 



about 3 /xG where it connects up with the disk of the galaxy. 



To get an idea of how the uncertainties in our parameter estimates affect our estimated 



12 Here a represents the spectral index of the non-thermal emission. 
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magnetic field strength, we varied each parameter individually by the estimated error in 
each parameter, calculated the difference between the new magnetic field estimate and the 
nominal estimate shown in Figure [9j and finally added the differences from each parameter 
variation in quadrature. We varied a by -0.3, Kq to 40 (a value typical for strong shocks), 
the distance by 10%, the diameter by 20%, the position angle by 3°, the center by 6" in each 
direction, and the inclination by 3°. The resulting calculation showed that the uncertainties 
on the total magnetic field strength in the central region and the outer region of the galaxy 
were 10 /zG and 3-4 /zG, respectively. This error is almost entirely due to the error in the 
random field. The error in the determination of the uniform field is only 0.5-2 /zG . If the 



Burn (1966) formula for q is used instead of the Sokoloff et al. (1998) formula for q, the 
uniform magnetic field is larger by a factor of 1.2 in the center of the galaxy and remains 
roughly the same at the edge and the random field remains the same everywhere. 

We also investigated how systematic uncertainties in the estimation of the synchrotron 
flux due to internal extinction in the central region of the galaxy will affect the estimated 
magnetic field. Assuming that the maximum internal extinction of 1.74 holds over the entire 
galaxy, we get a synchrotron flux down by a factor 0.6 near the main radio continuum peak, 
0.5 near the secondary radio continuum peak, and 0.2 in the region between the two peaks. 
Since the radio continuum emission is entirely random in the central region, 

B t oc Il ! ^- a) (22) 

and using the value of a determined above (-0.97), the maximum factor by which the mag- 
netic field strength will be decreased is 0.88 for the main radio continuum peak, 0.84 for the 
secondary radio continuum beam, and 0.67 for the region between the two peaks. 



4.3. Magnetic Field Structure 

The structure of NGC 1569's magnetic field reveals important clues about the interaction 
of the magnetic field with the other components of the interstellar medium and about the 
mechanisms at work in the galaxy that may have amplified the field. From our data set, 



we have two complementary measures of the structure of the magnetic field. In Figure 10 
the 3 cm polarization vectors rotated by 90° give an estimate of the orientation, but not 
direction, of the magnetic field in the plane of the sky. There is very little Faraday rotation 
at 3 cm, so the estimate of the magnetic field orientation at 3 cm should be accurate. The 
difference in the position angle of the polarized emission at different wavelengths gives the 
rotation measure 

= RM A 2 (23) 
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which depends on the thermal electron density and the magnetic field along the line of sight 



RM = 0.812 / B\\ n e ds 



(24) 



where RM is in rad m 2 , n e is in cm 3 , Bn is in fiG , and ds is in pc (Gardner & Whiteoak 



1966). The sign of the rotation measure gives the direction, not just the orientation of Bu. 
Positive rotation measures point towards the observer and negative away. In Figure 10 we 



plot the rotation measures between the 6 cm and 3 cm data with the estimated magnetic field 
vectors from the 3 cm data overlaid. We do not show the rotation measures between between 
the 20 cm and 13 cm data because there are only two patches of significant rotation measures. 
These patches are the size of the beam and are located at 4 h 30 m 41!41, 64 52'12"98 (-36 
rad m -2 ) and 4 h 30 m 45!17, 65°51'52 / . / (42 rad m -2 ). The rotation measures were determined 
using the MPS task COMB. 

Generally, polarization observations at two frequencies do not have enough information 



to accurately determine rotation measures (Reuter & Klein 1996). However, examining 



Figures [5] and [6j we find that the changes in angle within each pair of frequencies are small 
enough that we are able to determine the rotation measure between the frequencies in each 
pair with reasonable accuracy. For the polarization vectors to have rotated at least 180° 
between 6 cm and 3 cm, the rotation measure would have to be greater than 1230 rad m -2 , 
which is much larger than the rotation measures seen in other starbursting galaxies (see 



£4.4.3). Large changes in the distribution of the polarized emission between 3 cm and 20 cm 



that prevent the determination of rotation measures using at least three frequencies are a 



common feature of magnetic field observations of edge-on starbursting galaxies (Reuter et al. 



1994 Beck et al. 1994; Dahlem et al. 1997 Golla & Hummel 1994). 



From Figure [T0| one can see from the magnetic field vectors that, in general, the mag- 
netic field in the plane of the sky is roughly perpendicular to the disk of the galaxy. The 
magnitude of the rotation measures in the northeastern half of the galaxy are larger than 
those in the southwestern half of the galaxy. The relative magnitude of the rotation mea- 
sures between these two regions supports the geometry derived from X-ray measurements: 
the southern half of NGC 1569 is inclined towards us, while the northern half is inclined away 
( |Martin et al.|2002[ ). Therefore, the emission from the region north of the disk passes through 
a larger column of ionized gas and thus has a correspondingly larger rotation measure. In 
the southeastern portion of the galaxy, there are more rotation measures close to zero which 
means that, in general, the field we are seeing is mostly in the the plane of the sky and that 
the field along the line of sight is small. This overall magnetic field configuration means that 
the field can help transfer energy and momentum from cosmic rays to the thermal plasma 



(Breitschwerdt et al. 1991, 1993 Everett et al. 2008) helping to sustain the outflow of gas 



from NGC 1569. 
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The magnetic field structure of NGC 1569 is closely tied to the bubbles of ionized gas 
seen in Ha (see Figure W\ for a finding chart of the Ha bubbles identified in NGC 1569). 



In Figure [10] in the northeastern portion of the galaxy, we see a pair of oppositely directed 
rotation measures. The field along the line of sight in the easternmost region is directed 
toward the observer (yellow/green region in the northeastern part), while the field to the 
northwest of it is directed away from the observer (the blue region in the northeastern part). 
These rotation measures coincide with Ha bubble E. The combination of the orientation 
of the magnetic field vectors in the plane of the sky with the observed rotation measures 
suggests a magnetic loop being drawn out of the disk. The magnetic field vectors at 20 cm 
(see Figure [5] for the polarization vectors, which can be rotated by 90° to get an estimate 
of the magnetic field vectors) and the two small regions of significant rotation measures 
between the 20 cm data and the 13 cm mentioned earlier are associated with Ha Bubble B. 
The sign of the rotation measures between the 20 cm and 13 cm data is oppositely directed 
in the two regions, suggesting another field loop especially given the estimated orientation 
of the magnetic field in the plane of the sky. 

Looking closely at the rotation measures between 6 cm and 3 cm in the southern portion 
of the galaxy, one can see indications of other possible magnetic field loops associated with 
Ha bubbles A and F. The southern patch of rotation measures at 4 h 30 m 50 s shades from blue 
to green, which means that the field along the line of sight is going from being pointed away 
from the observer to pointed toward the observer. This region is associated with Bubble F. 
The patch of rotation measures at 4 h 30 m 42 s shows the opposite trend. The field goes from 
being almost completely in the plane of the sky (green) to pointed away from the observer 
(blue). The western arm goes from yellow where it connects to the disk of the galaxy (field 
pointed towards observer) to green (field primarily in the plane of the sky) to blue-green 
(field pointed away from observer). These regions are both associated with Bubble A. 

The bases of all the bubbles are around 30" which at the distance of NGC 1569 is a 
physical size of about 500 pc. This size is roughly double the size of the based of an average 



superbubble in the Milky Way (200 pc; Ferriere 2001 ) 



4.4. Implications of Derived Magnetic Field Strength and Structure 

44.1. Role of the Field m NGC 1569 

From our magnetic field estimates, we can calculate a typical pressure for the central 
region and halo of NGC 1569. Adopting 38 /iG for the central field and 12 /zG for the halo 
field (see Figure [9j, we get pressures (P/k) of 4.2 x 10 5 K cm -3 and 0.41 x 10 5 K cm -3 
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for the central magnetic pressure and the halo magnetic pressure. In Table [7| we compare 
these values with the pressures determined for other components of NGC 1569's ISM and 
the overall gravitational pressure. Where necessary, we have corrected the input values for 
the updated distance to NGC 1569. Except for the pressure of the diffuse ionized medium 
(DIG), we have not included any corrections for the highly uncertain filling factors of various 
component sj 1 ^ Note that we associate the broad component of the turbulent pressure with 
the DIG and the narrow component with the HII regions. Table [7] shows that the magnetic 
pressure in NGC 1569 is generally the same order of magnitude as most of the other sources 
of pressure in its ISM, but is also generally less than the other pressures in similar regions. 
The dominant pressure source in NGC 1569 is actually the turbulent pressure of the HII 
shells, which is not surprising given the level of disruption in this galaxy. The central x-ray 
pressure is almost as large as the gravitational pressure which explains the extended halo 
of hot gas around NGC 1569. The pressure balance in NGC 1569 is very different from 
the pressure balances in NGC 6946 QBeck||2007[ ) and M33 ( |Tabatabaei et al.||2008"| ). For 
these galaxies, the energy densities (i.e. pressures) of the magnetic field and the turbulent 
component are roughly comparable, while the pressure for the thermal gas is down by a 
factor of approximately 5. NGC 6946 and M33, however, are very different galaxies from 
NGC 1569. They are relatively quiescent, star-forming galaxies, unlike NGC 1569, which is 
in final throes of a starburst that has significantly affected its ISM. 



4-4- 2- Generating the Field 



Ferriere & Schmitt (2000, and references therein) and Gressel et al. (2008) have devel 



oped models in which galactic rotation provides the shear and supernova explosions drive the 
turbulence needed to generate a large scale magnetic field. In these models, one would expect 
the magnetic field to be correlated with the supernovae driven bubbles as seen in NGC 1569 
(see § 4.2). Note that this is a necessary, but not sufficient condition for a supernova driven 



dynamo to work. The calculations of Ferriere et al. and Gressel et al. are based on mean field 
a — bj dynamo models, which have been criticized on theoretical grounds (Kulsrud & Zweibel 



2008 Cattaneo & Hughes 2009). Nevertheless, whether or not mean field theory is valid, it 



seems likely that shear and turbulence are key ingredients in a galactic dynamo. Therefore, 
we estimate the strength of these two components (turbulence via supernova explosions and 
shear via differential rotation) for NGC 1569. 



13 Decreasing the filling factors from one, which is what we assume here, will increase the pressure by a 
factor of \jyfj . 
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First, we compare the supernova rate in the Milky Way to that of NGC 1569. The 
supernova rate from type II supernovae in the Milky Way is 27 Myr -1 kpc -2 (Ferriere 



200l[ ) J^| |Angeretti et al. (2005) find that the recent star formation rate for NGC 1569 is 
6.4 M yr" 1 kpc -2 . This star formation rate was determined using a distance estimate of 
2.9 Mpc. Increasing the distance used in the star formation rate estimate to the best current 
estimate (3.36 Mpc) will increase the rate by a factor of two to three and cause the latest 



epoch of star formation to end recently (Grocholski et al. 2008). Assuming a Salpeter IMF 



with 0.8 and 150 M as the minimum and maximum stellar masses and normalizing it to 
the total number of stars formed per year and per kpc 2 , we find that that the number of 
stars formed with masses greater than 8 M for NGC 1569 is 1.1 x 10 5 Myr -1 kpc -2 . With 
a duration of the most recent burst of star formation of 19 Myr ( Angeretti et al.|2005 ) and a 
typical lifetime for an 8 M star of 30 Myr, we get a supernova rate of 7 x 10 4 Myr -1 kpc -2 . 
This rate is a factor of 2500 higher than the supernova rate in the Milky Way. 



The shear is defined as 



d~R 



(25) 



where Q is the angular speed (v/R) and R is the radius. Given that the rotation velocity in 
the outer regions of galaxy disks is roughly constant, the shear can be estimated by 



d(v/R) 



dR 



R 2 ' 



(26) 



In the Milky Way, the angular speed is approximately 220 km s" 1 at a radius of 8.0 kpc, so 
the shear is 3.4 km s -1 kpc -2 . In NGC 1569, the galaxy has a rotation speed (corrected for 
an inclination of 60°) of 46 km s -1 at a radius of 1.4 kpc (Miihle et al. 2005), which gives a 
shear of 33 kms -1 kpc -2 , which is a factor of 10 larger than the shear in the Milky Way. 

Since both the shear and the supernovae rate are much greater than in the Milky Way, 
the fidicual model for this type of dynamo, it is quite possible for a supernova-driven dynamo 



like those described in Ferriere & Schmitt (2000) and Gressel et al. (2008) to be at work in 



NGC 1569. Detailed modeling, however, is needed to determine whether or not this type of 
dynamo is at work in NGC 1569. Note that the supernova-driven dynamo described above 



is different that the related mechanism modeled described by Parker (1992) and modeled 



by Otmianowska-Mazur et al. (2009). In that mechanism, the magnetic buoyancy driven by 



cosmic rays lofts the field loops into the halo of the galaxy. 



including type I supernovae increases the total rate by 15-20%. 



4-4-3- Comparison to Magnetic Fields of Other Starbursts 



Table [8] summarizes the estimated magnetic field strengths in other starburst galaxies. 
Overall, we find that our estimate for the magnetic field strength of NGC 1569 is in the 
middle of the distribution of magnetic field strengths for starbursting galaxies with NGC 
1569's field being greater than that of NGC 253 and NGC 4666, but less than that of M82 
and NGC 1808. It is important to note, however, that these galaxies are all much larger than 
NGC 1569, so for the field in NGC 1569 to be even comparable to that in larger galaxies has 



interesting implications for dynamo theory. Thompson et al. (2006) argue that the magnetic 



field strengths of starbursts like NGC 253, M82, and NGC 1808 are underestimated by the 
minimum energy magnetic field estimates based on synchrotron emission and argue for much 
greater fields strengths in starbursting galaxies based on equipartition of the gravitational 



and magnetic pressures. However, observations of Zeeman splitting in (U)LIRGS by Ro- 



bishaw et al. (2008), which yield direct measurements of the magnetic field strength, agree 



roughly with the values derived from synchrotron emission. 

Comparing the magnetic field structure we derive for NGC 1569 (see Figure [9]) with the 
magnetic fields structures derived for other starburst galaxies (NGC 253, M82, NGC 1808, 
NGC 4569, and NGC 4446), we can see several important differences. First, the polarization 
features in NGC 1569 are patchier than in any of the other galaxies. There are two probable 
causes for this patchiness: 1) the magnetic field of NGC 1569 is much more turbulent and 2) 
the depolarization effects in NGC 1569 are much more severe. Given the turbulent nature of 



the Ha in the central regions of NGC 1569 (Westmoquette et al. 2008), it is likely that the 



field is very turbulent in this region as well, although the large electron density here makes 
it possible for Faraday depolarization to mask any uniform field component. However, in 
the outer regions of this galaxy depolarization effects due to the Ha filaments are likely to 
dominate over turbulence given the organization imposed on the ISM by the gas outflow. In 
addition, if the depolarization in the outer regions was due to turbulence it would mean that 
the turbulence would have to vary in the halo of the galaxy, which is at odds with the result 



of Westmoquette et al. (2008), who found that the warm ionized gas in the outer regions 
has more or less constant velocity width between 20 to 70 km s _1 . The second difference 
between the previous observations and our observations of NGC 1569 is that the ratio of the 
length of the polarized arm to the major axis of the radio continuum emission associated 
with the galactic disk is very high (~ 0.5). Of the previous observations of starbursting 



galaxies, only NGC 4569 (Chyzy et al. 2006), whose interstellar medium has been shaped 



by its location in the Virgo cluster, has such a large ratio of radio continuum feature length 
to length of the major axis of the radio continuum emission associated with the disk. 



Rotation measures have been determined for M82 (Reuter et al. 1994), NGC 253 (Beck 



etlHl994fr , and NGC 4446 ( |Dahlem et al|l997fr . Comparing our derived rotation measures 



in Figure 10 to those derived between 6 cm and 3 cm for M82 and NGC 253, we see that the 
average rotation measure between these two frequencies is more or less the same for all three 
galaxies (~ ±200 radm -2 ). The rotation measures for NGC 4666 were determined between 
20 cm and 6 cm and their magnitudes correspond well to the rotation measures determined 
here between 20 cm and 13 cm (~ ±40 radm -2 ). 



5. A Multiwavelength Picture of the Interstellar Medium of NGC 1569 



NGC 1569 is a well-studied galaxy. In this section, we compare our observations with 
those at other wavelengths to examine the larger picture of the complex interstellar medium 
of NGC 1569. We start with the Ha emission from NGC 1569, which is characterized by 
long Ha filaments extending along the minor axis of the galaxy. The filaments were first 



discussed extensively by 


Hodge 


(1 


974) 


and were the subject of much subsequent work (e.g. 


de Vaucouleurs et al. |1974 Waller 


199] 


L Hunter et al. |1993[ |Heckman et al. 1995 Hunter & 


Gallagher||1997 


)■ 


Westmoquette et al. ( 


2007alb 


) show that in the central region the emission 



lines have two components: a narrow-line component (FWHM < 70 km s _1 ) that they 
associate with the turbulent ISM and a broad-line component (FWHM ~ 150 km s -1 ) that 
they associate with emission from a turbulent mixing layer between cool gas clumps and the 
hot gas from surrounding star clusters. The broad-line component of the emission is confined 



to the central region of the galaxy (Westmoquette et al. 2008). The outer regions of NGC 



1569 show clear signs of an outflow of gas from the center (Westmoquette et al. 2008). This 
result agrees with previous results from long-slit spectra ( Tomita et al.||1994 Heckman et al. 
1995l|Martin||1998l). 



Figures [5] and [6] compare the distribution of the Ha emission to the radio continuum 
emission and were partially discussed in § |3.2| and |4.3| Here we discuss three additional 



points. First, the radio continuum arm becomes increasingly aligned with the Ha arm as 
the wavelength decreases. This fits well with the region on the outside of the Ha arm 
being synchrotron dominated, which is what one would expect from our thermal fraction 
estimates in § 44 Second, the uniform component of the magnetic field in this region may 
be at least partially due to compression by the expanding Ha bubble in that region. In 
Figure |9j the region of uniform magnetic field lies outside of (to the west) of the Ha arm. 
Finally, the presence of a magnetic field in the Ha arm region may explain why the western 
Ha arm is thicker than the other Ha features in this galaxy (see Figures [5] and [6]). The 
field here could provide an additional source of pressure, working against any compression 
forces, thereby increasing the width of the filament compared to the filaments elsewhere in 



-29- 



the galaxy. However, the effect of the magnetic field on the arm is depends on a variety of 
factors including the magnetic field strength in the arm, the variation of thermal pressure 
with height above the galaxy disk, and the effect of the starburst on the distribution of ISM 
of NGC 1569. In addition, our data only gives the integrated magnetic field - we cannot 
determine what the field is in the Ha arm, only the integrated field in that region, which 
is dominated by depolarized effects. Detailed modeling is the only way to disentangle these 
effects. 

X-ray observations of NGC 1569 reveal that the galaxy is surrounded by a halo of hot 
gas bounded by Ha filaments (Heckman et al. 1995 della Ceca et al. 1996). High spatial 



resolution observations by Martin et al. (2002) have shown that the bright X-ray emission 



near the filaments is most likely due to shocks generated by the outflow of gas from NGC 
1569. The measured ratio of a elements to iron is between 2.5 and 4.8 times the solar 



ratio, indicating that the wind is laden with the products of type II supernovae (Martin 



et al. 2002). We compare our radio continuum observations with the adaptively smoothed 
X-ray map of Ott et al. (2005a) in Figure 11 To compare the extents of the X-ray and 



radio continuum maps, we have derived scale lengths for the radio continuum emission by 
azimuthally averaging the emission as a function of radius from the center of the emission 
using IRING in AIPS. We find that the scale lengths are 376.1 pc (20cm), 310.9 pc (13 cm), 



236.5 pc (6cm), and 258.7 pc (3cm). Ott et al. (2005b) found that the scale length for the 
X-ray emission is 208 pc, so the radio continuum emission extends further than the X-ray 
emission at all wavelengths. We suggest that the magnetic field in the radio continuum arm 
may confine the hot gas on the western edge of the galaxy and funnel the gas toward the 
halo of the galaxy. Similar effects are seen in magnetohydrodynamic simulations of the ISM 



by de Avillez & Breitschwerdt (2005). 



The neutral hydrogen distribution of NGC 1569 is quite complex (cf. Stil & Israel 



2002 


Miihle et al. 


2005 


and 


Miihle et al. 


2006) 



The HI distribution has an arm to the 
west of the galaxy, but this feature is distinct from the arm seen in X-ray, Ha, and radio 



continuum. Miihle et al. (2005) suggest that the HI arm may be the remnant of a primordial 
gas cloud that interacted with NGC 1569 and started the most recent burst of star formation 
in this galaxy; the HI arm is likely still in the process of accreting onto the disk of the 



galaxy. In Figure 12, we compare the neutral hydrogen distribution to our radio continuum 
observations. The location of the polarized western arm south of the HI arm also suggests 
an interaction between the infalling HI arm and the galaxy itself. We speculate that the 
strong magnetic field associated with the Ha arm is the result of the collision of the HI arm 
with the expanding bubble (Bubble A in Westmoquette et al.||2008 and Filament 6 in Hunter 



et al. 1993). The western Ha arm is morphologically very different from the other filaments 



in the galaxy and given the confluence of interesting Ha, HI, and radio continuum features 
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in this region, it would be surprising if they were not related. 



The extended radio halo around NGC 1569 may not be unique to this galaxy. In general, 
dwarf galaxies have very shallow potential wells, so any significant star formation episodes 
will tend out drive material out of their potential wells. This mechanism has been postulated 
as a way to enrich the intergalactic medium ( Mac Low fc Ferrara||1999 ), but it could also be 



a mechanism to magnetize the intergalactic medium (Kronberg et al. 1999 Volk & Atoyan 



2000} |Bertone et al.|2006| >. If these halos are like the radio continuum halo of NGC 1569, then 
they will be dominated by synchrotron emission. Therefore, the best wavelength regime to 



search for them is the low frequency radio regime with telescopes like LOFAR (Klein 2008). 



The detection of such halos in galaxies other than NGC 1569 would provide important clues 
about the magnetization of the universe and the escape of cosmic rays from the halos of 
galaxies. 



6. Summary and Conclusions 



NGC 1569 provides a crucial test of the importance of magnetic fields in starbursting 
dwarf galaxies. We present the first observations to date of the polarized radio continuum 
emission at 20 cm, 13 cm, 6 cm, and 3 cm from the post-starburst dwarf irregular NGC 1569. 
We detect radio continuum emission from NGC 1569 at all four wavelengths. Our observa- 
tions are in agreement with previous high-resolution radio continuum observations of NGC 



1569 (Seaquist & Bignell 1976 Condon 1983 Klein & Grave 1986; Israel & de Bruyn 1988 



Lisenfeld et al.|2~004 ). We confirm the discovery of an extended radio continuum halo, which 



was first imaged by Condon (1983) and discussed by Israel & de Bruyn (1988), and for the 



first time clearly image the western arm at wavelengths shorter than 20 cm. In general, the 
spectral indices derived for this galaxy get steeper as one moves along the minor axis of 
the galaxy away from the disk. However, there are two filamentary regions of flat spectral 
indices extending to the edge of the map. The spectral index trends shown here agree with 



those found by Lisenfeld et al. (2004) supporting their theory that there is a convective wind 



at work in the halo of this galaxy. 

We detect strong polarized emission at all four frequencies. The polarized emission 
at 6 and 3 cm is spread over most of the galaxy, while the polarized emission at 20 cm is 
much spottier, and there are only a few small patches of polarized emission at 13 cm. The 
polarized emission at 20 cm, 6 cm and 3 cm is associated with regions identified as bubbles 
from Ha observations. There is little polarization in the central region of NGC 1569, where 
the ISM is very turbulent and the largest thermal electron column density can be found. 
The Ha filaments in the halo also depolarize the emission. 
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In order to determine the magnetic field strength, we estimate the thermal contribution 
to the radio continuum emission using a flux-calibrated Ha image. We find that the thermal 
contribution peaks at 40-50% in the center of the galaxy and falls off quickly. The halo is 
dominated by synchrotron emission. The total magnetic field strength in the central regions 
is 38 /zG and approximately 10-15 /iG in the halo. The field is largely random near the 
center; the uniform field is strongest in the halo. The magnetic field structure of NGC 1569 
is shaped by the outflow of gas from this galaxy drawing field loops out of the disk. We show 
that the magnetic field structure of NGC 1569 is closely aligned with the Ha bubbles. 

We find that the magnetic pressure is the same order of magnitude as other components 
of the pressure in the ISM, but is generally less than the pressures of other components in 
similar regions by a factor of a few. Therefore, the magnetic field does not dominate the 
ISM of NGC 1569, although it clearly plays an important role. 

The close correspondence between the Ha bubbles and the magnetic field structure 
suggests that the interstellar medium of NGC 1569 contains magnetic bubbles similar to 



those modelled by Ferriere et al. (1991). The high supernova rate and strong shear in the 



galaxy compared to the Milky Way suggests that these bubbles might combine to give a 



supernova/superbubble driven dynamo like that of Ferriere & Schmitt (2000) and Gressel 



et al.| ((2008 ) . However, detailed modeling is necessary to say definitely whether or not this 
type of dynamo is at work in NGC 1569. 

Finally, we compare our data to observations of the ISM of NGC 1569 at other wave- 
lengths. The turbulent Ha kinematics of the central portion of the galaxy are reflected in a 
primarily random field, while the outflow kinematics in the halo are reflected by a uniform 
field. The magnetic field in the Ha arm may be increasing the width of this feature by 
providing an additional source of pressure as well as channeling hot gas out of the plane. We 
suggest that the infall of the western HI arm into the disk of NGC 1569 and the outward 
expansion of the bubble associated with the western Ha arm compressed the ISM in this 
region and produced the strong uniform magnetic field seen here. In short, our observations 
reveal that the magnetic field of NGC 1569 is shaped by the outflow of gas from the galaxy 
and its interaction history and that the magnetic field may be playing an important role in 
directing the outflow of gas from the galaxy. 

Radio continuum halos generated by the outflow of gas from galaxy may be common in 
dwarf starburst galaxies. These halos provide a mechanism for magnetizing the intergalactic 
medium ( Kronberg et al.|1999 Volk fc Atoyan|2000 Bertone et al.||2006 ). Surveys of edge-on 
dwarf starbursts with future low-frequency radio telescopes like LOFAR will help constrain 
the number of dwarf starbursts with synchrotron-dominated radio halos and provide impor- 
tant clues about the magnetization of the universe and the escape of cosmic rays from galaxy 
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halos. 
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A. Deriving the Geometric Corrections 

We assume that the geometry of NGC 1569 can be approximated by an infinite cylinder 
with the vector normal to the galaxy disk inclined from the line of sight by % degrees (i = 0° 
is face-on, i = 90° is edge-on) and rotated on the sky by a position angle (measured from 
north and increasing counterclockwise, i.e, towards the east) PA g . The geometry for NGC 
1569 and the relevant variables are shown in Figures 13 and 14 Below we derive the relevant 
quantities needed for our calculation of the magnetic field strength. 



A.l. Deriving the Line of Sight Distance 

The line-of-sight length (I) through an inclined (infinite) cylinder is 

l = H# : W (A1) 

smz 

where d is the distance from the central axis of the cylinder, R is the radius of the cylinder, 
and i is the inclination from the line of sight. If the cylinder is rotated about its minor axis 
in the plane of the sky, where the angle between the RA axis and the disk of the galaxy 
(minor axis of cylinder) is (3, d becomes 

d = D-^— [(RA pix - RA ) cos (Dec ) cos (5 - (Dec pix - Dec ) sin (5] (A2) 
180 
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where D is the distance to the galaxy, (RA pix , Dec pix ) are the coordinates of a particular 
pixel, and (RA , Deco) are the coordinates of the center of the galaxy. The variable (3 is 
defined as positive for counterclockwise rotations and negative for clockwise rotations and is 
90 + PA g for PA g < and PA g - 90 for PA g > 0. If /3 is 0°, then d is just 

d = L>-^-(RA pix - RA ) cos (Dec ) (A3) 



A. 2. Deriving c 4 for our geometry. 

To deal with the position angle of the galaxy, we correct the observed polarization vector 
position angles (PAb,o)> which are on the usual position angle system where PA g = 0° is 
north, PAg = 90° is east, and PA g = —90° is west, so that they are centered on a coordinate 
system where a position angle of 0° is along the axis of the cylinder and a position angle of 
90° is pointing toward the left (eastward) and a position angle of —90° is pointing toward 
the right (westward). For a negative position angle, the correction is 

PA B = - (90° + PAg) + PA B ,o (A4) 

and for a positive position angle 

PA B = (90° - PAg) + PA B ,o (A5) 

We rotate any resulting PA B outside the range of [90,-90] by 180°. 

In the frame of the cylinder, there are two components to the observed magnetic field 
vector (B±): B x and B Z} ±. The B x components are not affected by the inclination of the 
galaxy because they are perpendicular to the inclination axis. In the z direction, however, 
we only see the field strength B z projected onto the plane of the sky (B z ^±). The observed 
magnetic field strength in the z direction is 

B Zt± =B z sim (A6) 

where B z is the magnetic field strength in the z direction and i is the inclination. Now the 
observed magnetic field strength (B±) is also equal to 

B z , ± = B ± cos (PA B ) (A7) 



Therefore, 



B z = B ± cos (PA B ) / sin i. 



(A8) 
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Finally, the total magnetic field (B t ) is 

(B ± sin(PA B )) 



2 , B± cos (PA B ) 
sini 



?.l 1/2 



(sin(PA B )) 2 + 



2 f cos (PAb^ 



sin? 



1/2 



Beck & Krause (2005) define C4 as 



C4 



(1-a) 



so our C4 is then 



c 4 



(sin (PA B )) 2 + 



2 / cos (PAb^ 



sinz 



-(l-a)/2 



for a completely regular field. 



(A9) 
(MO) 

(AH) 



(A12) 



(A13) 
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Table I. Summary of WSRT 20 cm Observations 





5-6 Nov 
2000 


6-7 Nov 
2000 


Array Configuration 


36m 


72m 


IF1 (GHz) 


1.345 


1.345 


IF2 (GHz) 


1.355 


1.355 


IF3 (GHz) 


1.365 


1.365 


IF4 (GHz) 


1.375 


1.375 


IF5 (GHz) 


1.385 


1.385 


IF6 (GHz) 


1.395 


1.395 


IF7 (GHz) 


1.405 


1.405 


IF8 (GHz) 


1.415 


1.415 


Bandwidth per IF (MHZ) 


10 


10 


Flux Calibrator 


0137+331 


1331+305 


Instrumental 


0137+331 


0542+498 


Polarization 






Absolute 


2253+161 


1331+305 


Polarization Angle 







Table 2. Final Images 





Frequency 




Beam 


PA 


LAS a 








Data Name 


GHz 


Weighting 


// 


O 


/ 


fijy beam -1 


fijy beam -1 


/ijy beam -1 


20cm 


1.38 


Natural 13 


12.88 x 14.19 


6.5 


21 


16.5 


12.6 


13.2 


13 cm, robust=0 


2.27 


Robust = 


10.19 x 9.49 


-9.46 


13 


19.6 


17.2 


17.3 


13 cm, 13" 


2.27 


uvrange=(0.165,13) 


13.44 x 12.96 


-9.60 


13 


24.8 


19.4 


19.3 


6 cm, robust=0 


4.86 


Robust = 


4.42 x 4.08 


89.77 


c 


8.5 


7.6 


7.6 


6 cm, 13" 


4.86 


uvrange=(0.165,13) 


12.64 x 11.66 


-88.20 


c 


14.0 


10.8 


10.7 


3 cm, robust=0 


8.46 


Robust = 


8.85 x 6.82 


-78.31 


c 


11.1 


8.0 


8.0 


3 cm, 13" 


8.46 


uvrange=(0.165,13) 


12.65 x 11.26 


-86.44 


c 


17.6 


10.8 


11.0 



a Largcst Angular Scale 

b All fixed-moveable baselines selected. 

c This image is the result of combining VLA and Effclsbcrg observations, so it includes all spatial frequencies greater than the frequency of 
the synthesized beam. 
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Table 3. Summary of WSRT 13 cm Observations 





4-5 Dec 


18-19 Jan 




2002 


2003 


Array Configuration 


72m 


36m 


IF1 (GHz) 


2.200000 


2.20000 


IF2 (GHz) 


2.218125 


2.218125 


IF3 (GHz) 


2.236250 


2.236250 


lr 4 (Cjrlz) 


2.254375 


2.254375 


IF5 (GHz) 


2.272500 


2.272500 


IF6 (GHz) 


2.290625 


2.290625 


IF7 (GHz) 


2.308750 


2.308750 


IF8 (GHZ) 


2.326875 


2.326875 


Bandwidth per IF (MHz) 


20 


20 


Bandpass Calibrator 


1331+305 


1331+305 


Flux Calibrator 


1331+305 


1331+305 


Instrumental 


0834+555 


0834+555 


Polarization 






Absolute 


1331+305 


1331+305 


Polarization Angle 






R-L Angle 


66 


66 



Table 4. Summary of VLA 6 cm and 3 cm Observations 





11 Aug. 


13 Aug. 


17 Aug. 


27 Aug. 


19 Sept. 


23 Sept. 


9 Aug. 


10 Aug. 




2000 


2000 


2000 


2000 


2000 


2000 


2001 


2001 


Program Number 


AM643 


AM643 


AM643 


AM643 


AM643 


AM643 


AM694 


AM694 


Array 


D 


D 


D 


D 


D 


D 


C 


C 


IF1 (GHz) 


8.4351 


4.8851 


8.4351 


8.4351 


8.4351 


8.4351 


4.8851 


4.8851 


IF2 (GHz) 


8.4851 


4.8351 


8.4851 


8.4851 


8.4851 


8.4851 


4.8351 


4.8351 


Bandwidth (MHz) 


50 


50 


50 


50 


50 


50 


50 


50 


Flux Calibrator 


0521+166 


1331+305 


0521+166 


1331+305 


1331+305 


1331+305 


0137+331 


0137+331 






0521+166 






0521+166 


0521+166 






Phase Calibrator 


0244+624 


0614+607 


0244+624 


0614+607 


0244+624 


0244+624 


0614+607 


0614+607 












0614+607 


0614+607 


0217+738 


0217+738 


Instrumental 


0244+624 


0319+415 


0244+624 


0614+607 


0244+624 


0244+624 


0614+607 


0614+607 


Polarization 




0614+607 






0614+607 


0614+607 


0217+738 


0217+738 


Absolute 


0521+166 


1331+305 


0521+166 


1331+305 


1331+305 


0521+166 


0521+166 


1331+305 


Polarization Angle 


















R-L Angle 


-20° 


66° 


-20° 


66° 


66° 


-22° 


-24° 


66° 



Table 5. Summary of Effelsberg 100-m Observations 



Frequency 


Bandwidth 




Beam 


Map Size 


Observing Time 


Noise 


GHz 


MHz 


Number of Feeds 






hr 


mjy beam -1 


4.85 


500 


2 


2.7 x 2.5 


10 x 10 


5 


2.5 


10.45 


300 


4 


1.1 x 1.1 


12 x 8 


26 


0.5 
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Table 6. Parameters of Previous Observations 





Wavelength 


Frequency 


Resolution 


Largest Angular Scale 


Sensitivity 


Paper 


CHI 


GHz 


it 




mjy beam -1 


Seaquist & Bignell (1976) 


11.1 


2.7 


7.7 x 5.5 


1 


0.71 mjy beam 




3.7 


8.1 


3.0 x 2.1 


1.3 


0.23 mjy beam 


Condon (1983) 


21.4 


1.4 


~ 15 - 20 


15 


0.04 mjy beam -1 




21.4 


1.4 


5.9 x 5.8 


2 


0.06 mjy beam 


Klein & Grave (1986) 


1.2 


24.5 


29.7 


a 


1.7 mjy beam -1 


Israel & dc Bruyn (1988) 


50 


0.6 


56 x 62 


48 


b 




21.4 


1.4 


24 X 27 


14 


0.65 mjy beam -1 




6 


5.0 


6.9 x 7.6 


6 


b 


Lisenfeld et al. (2004) 


20 


1.5 


6 


2 


0.12 mjy beam 




6.1 


1.9 


6 


2 


0.13 mjy beam -1 




3.6 


8.4 


6 


2 


0.09 mjy beam -1 




1.9 


15.4 


6 


2 


0.15 mjy beam - 



a These are single-dish observations, which are sensitive to all spatial scales greater than their resolution. 
b Sensitivities not given. 



Table 7. Pressures of Various Components of the ISM of NGC 1569 



Pressure 



Component 


Input Values 


Equation 


Reference 


10 s K cm -3 


Magnetic Field (central) 


B 


= 38 imG 


B 2 /(8Tvk) 


1 


4.2 


Magnetic Field (halo) 


B 


= 12 aiG 


B 2 /(8TTk) 


1 


0.41 


Hot Gas (central) 


n e = 0.049 cm" 


" 3 , T = 7.23 x 10 6 K 


2n e T 


2 


7.2 


Hot Gas (halo) 


n e = 0.028 cm" 


~ 3 , T = 7.23 x 10 6 K 


2n e T 


2 


2.0 


HII Regions 


n e ~ 57 cm - 


3 , T = 1.0 x 10 4 K 


2n e T 


3 


11 


Diffuse Ionized Gas (DIG) 


equilibrium between hot gas and DIG 




4 


1-6.3 


Warm Gas (turbulent, broad) 


n e ~ 0.1 cm -3 


i v turb = 150 km S 1 


n e m H v 2 urb /(2k) 


5 


1.4 


Warm Gas (turbulent, narrow) 


n e ~ 57 cm -3 


i Vturb = 70 km s 


n e m HV 2 urb /(2k) 


3, 5 


170 


Gravity 


Vrad = 65.5 km s" 


_1 , i = 63°, r = 1.4 kpc 


3(iw/ sin i) 4 /(20Trr 2 Gk) 


6 


8.1 



References. — (1): This work, (2): |Ott et al.| ( |2005b| l, (3): |Martin| ( [l997| l , (4): |Martin| ( [1999^ , (5): |Westmoquette et al.| | |2008| l, (6): 
|Muhle et al] (2005) 
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Table 8. Magnetic Field Strengths in Starburst Galaxies 





B t 


Bt, center 


Bt,halo 




Galaxy 








Reference 


M82 


50 






1 


NGC 253 




17 


7 


2 


NGC 4666 




14.4 


7.1 


3 


NGC 1808 




75 


10-35 


1 


NGC 1569 




38 


10-15 


5 



a These are estimates of the magnetic field strength based 
on equipartition between gravitational and magnetic pres- 
sures from |Thompson et al.| ( |2006[ > . 



Referenc es. — (1): |Klein et al.| l |1988|l, (2): |Beck et al.| 
jl994| , (3): |Dahlem et al.| ( |1997^ (4): pahlem et a.1. | ( | 1990| > , 



(5): This work 
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Fig. 1. — 20 cm, 13 cm, 6 cm, and 3 cm total intensity, robust=0 weighted images. The 
contours start at 3a and increase by factors of 2 (see Table [2] for the noise levels for each 
image). The beam for each observation appears in the lower left corner. See Figure [2] for a 
more detailed view of the 6 cm and 3 cm emission. 
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Fig. 2. — Zoomed-in version of the 6 cm and 3 cm total intensity, robust=0 weighted images. 
The first contour is at 3cr and subsequent contours increase by factors of two (see Table [2] 
for the noise levels for various images). The beam appears in the lower left corner of each 
panel. 



- 46 - 



Wavelength (cm) 
1000 100 10 1 




100 1000 10000 

Frequency [MHz) 



Fig. 3. — Integrated flux density spectrum for NGC 1569. The small dots are from the 
compilation of Lisenfeld et al. (2004) and open circles are the total flux in our final images. 
The error in the total flux of our final images is less than the size of the symbol. The 
open triangles are the total flux densities determined from our Effelsberg observations. The 
8 GHz flux was determined by scaling our 10 GHz Effelsberg observations to 8 GHz assuming 
a spectral index of 0.47. The error bars show la errors. The line shows a power law with a 
slope of a = 0.47. 
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Fig. 4. — Spectral index map between 6 cm and 3 cm data. For this map, both data sets 
were imaged with the same u-v range (0.74 to 20 kA) and smoothed to the same beam size 
(9.34" by 7.54", shown in the lower left corner). Regions where the formal error on the 
spectral indices exceeded 0.2 were blanked. The contours show the 3 cm, robust=0 image 
and start at the 3<r noise level and increase by a factor of two for each subsequent contour. 
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Fig. 5. — E- vectors showing the direction and distribution of the polarized intensity at 
20 cm and 13 cm with contours showing the total intensity radio continuum emission at 
each wavelength (robust=0 weighting) and a grayscale image showing the distribution of 
the Ha emission. A polarization vector with a length of 5" has a polarized intensity of 51.1 
/iJy beam -1 . The contour levels are for the robust=0 image in each case. They start at 
3a and increase by a factor of 2 for subsequent contour (see Table [2] for the noise levels for 
each image). The central portion of NGC 1569 has been saturated in the grayscale image 
to better show the faint, extended emission. The beams are indicated the lower left hand 
corner of each panel. The line in the greyscale image roughly running from RA 04 h 30 m 45 s 
to 04 h 30 m 22 s at a declination of 64°47' is due to the edge of the CCD. 
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Fig. 6. — E- vectors showing the direction and distribution of polarized intensity at 3 cm and 
6 cm with contours showing the total intensity radio continuum emission at each wavelength 
and a greyscale image showing the distribution of the Ha emission. We show the 6 cm and 
3 cm data with 13" resolution because these were the data with the most polarized intensity. 
The contours start at 3<r and increase by a factor of 2 for each subsequent contour. A 
polarization vector with a length of 1" has as polarized intensity of 12.8 /xJy beam -1 . The 
beams are indicated in the lower left corner of each panel. 
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Fig. 7. — An Ha image image from Hunter & Elmegreen (2004) with outlines indicating the 
Ha bubbles from Martin (1998) and Westmoquette et al. (2008) and contours showing the 
3 cm radio continuum emission with 13" resolution. The contours start at 3a and subsequent 
contours increase by factors of two. See Table [2] for the noise levels in each image and Figure 
2 in Westmoquette et al. (2008) for the correspondence between this nomenclature and the 
nomenclature of Hunter et al. (1993). 
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Fig. 8. — Fraction of the total emission at 20 cm, 13 cm, 6 cm, and 3 cm that is thermal. 



The thermal emission was estimated using an Ha image of the galaxy (see § 4.1). Overlaid 



are the contours of total intensity at each wavelength with 13" resolution. The contours 
start at 3cr and increase by a factor of 2 for each subsequent contour (see Table [2] for the 
noise levels for each of the images). The beam is shown in the lower left hand corner of 
each panel. The thermal fraction image at 20 cm does not extend to the boundaries of the 
20cm emission because the 20 cm radio continuum emission extends beyond the Ha image 
boundaries, particularly on the east and north sides of the galaxy. 
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Fig. 9. — Magnetic field strength of NGC 1569 using the 3 cm data. The color plots are, 
clockwise from upper left, total magnetic field strength (in fiG), uniform magnetic field 
strength (in fiG), q, which is the ratio of the random magnetic field strength to the uniform 
magnetic field strength in the plane of the sky (this ratio is usually greater than one since 
the random field is usually much stronger than the uniform field) , and random magnetic field 
strength (in fiG). See [4.2 for the assumptions involved in the derivation of these quantities. 
The total, uniform, and random magnetic field strength have been corrected for the assumed 
geometry of NGC 1569 (see Figures 13 and 14). The 3 cm total intensity emission with 13" 
resolution is shown by contours. The contours are 3, 6, 12, 24, 48, 96, 192, 384, 768, and 
1536 times the la noise level. 
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Fig. 10.— Rotation measures between the 14" resolution 6cm and 3cm data (color), 
magnetic field vectors derived from the 3 cm data, and the contours of the 14" resolution 
3cm data overlaid on Ha image courtesy of D. Hunter (Hunter & Elmegreen 2004). The 
color bar is in units of 10 3 rad m~ 2 , a 1" polarization vector has a polarized intensity of 
12.8 /zJy beam -1 , and the contours are 3, 6, 12, 24, 48, 96, 192, 384, 768, and 1536 times 
the la noise. 
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Fig. 11. — Adaptively smoothed, total (0.3 keV < E < 8.0 keV) X-ray emission image from 



Ott et al. (2005a) with 20 cm, 13 cm, 6 cm, and 3 cm contours overlaid. The radio continuum 



images all have beams with approximately the same size as the 20 cm beam. The contours 
start at 3 times the la noise level in each band and each subsequent contour increases by 
a factor of 2. The beam for each wavelength is in the lower left of each panel. The region 
of dark gray in the lower left hand corner of each panel is not emission from NGC 1569 but 
instead is due to edge effects from the adaptive smoothing process. 
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Fig. 12. — Neutral hydrogen emission from Miihle et al. (2005 , 2006 ) with 20 cm, 13 cm, 6 cm, 
and 3 cm contours overlaid. The radio continuum images all have beams with approximately 
the same size as the 20 cm beam. The contours start at 3 times the la noise level and 
increase by a factor of two for each subsequent contour. The beam for the radio continuum 
images is shown the lower left of each panel. 
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Fig. 13. — Geometry assumed for NGC 1569 in plane of sky. 




Fig. 14- 



Geometry assumed for NGC 1569 along line of sight 



